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Abstract 

Sand-moving winds blowing from a constant direction in an area of high sand availability form transverse dunes, which 
have a fixed profile in the direction orthogonal to the wind. Here we show, by means of a linear stability analysis, that 
transverse dunes are intrinsically unstable. Any along-axis perturbation on a transverse dune amplify in the course of 
dune migration due to the combined effect of two main factors, namely: the lateral transport through avalanches along 
the dune's slip-face, and the scaling of dune migration velocity with the inverse of the dune height. Our calculations 
provide a quantitative explanation for recent observations from experiments and numerical simulations, which showed 
that transverse dunes moving on the bedrock cannot exist in a stable form and decay into a chain of crescent-shaped 
barchans. 
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1. Introduction 

Sand dunes are widespread on Earth deserts and coasts, 
and are also found on Mars, Venus, Titan and even on the 
h bottom of rivers [H, 0, H, B H S 0| ■ Since the pioneering 
Q work by tlL insi ghts from field and experimental works 
, Q;|, |, [lHjlfl^^ well as numerical simu- 
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our knowledge on the physics of sand transport and dune 
formation. 

Dunes form wherever sand is exposed to a wind that is 
strong enough to put grains into saltation — which con- 
sists of grains jumping in nearly ballistic trajectories and 
ejecting new par ticles upon collision with the sand bed 
El, H, IM lH S|- The main factors controlling dune 
morphology are the wind directionality and the amount 
of loose sand available for transport [27|. While elongat- 
ing, longitudinal seif dunes and accumulating star dunes 
develop under bi- and multidirectional wind regimes, re- 
spectively [8], in areas where the wind direction is nearly 
constant, two types of migrating dune may occur: 

• barchan dunes, which have a crescent shape and two 
limbs pointing in the migration direction (Fig. 
they occur in areas where the amount of sand is not 
sufficient to cover the ground [2 

0; 

• transverse dunes, which propagate with nearly in- 
variant profile orthogonally to a fixed wind direction 



Email addresses: hygorOfisica.ufc.br (Hygor P. M. Melo), 
parteli@fisica.ufc.br (Eric J. R. Parteli), soares@fisica.ufc.br 
(Jose S. Andrade Jr.), hansSif b .baug. ethz . ch (Hans J. Herrmann) 



(Fig. [TJd) ; they form when the amount of sand on the 
ground is high 35 , 36 , 37 , 3|| . The transverse dune is 
the simplest and best-understood type of dune, and 
occurs on all planetary bodies where dunes have been 
detected 0- 

As a matter of fact, barchans are observed in bedrock 
areas of Earth's desert and coastal environments [3j| or 
on the floor of Martian craters [40( where sand availability 
is low. Water tank experiments of dune formation under 
unidirectional stream also produce barchan dunes if the 
sand cover is incipient or there is no input of sand 41, 

t. These experiments, as well as numerical simulations 
, 43, suggest that transverse dunes cannot exist 
in a stable form on the bedrock: a transverse sand ridge 
evolving on the bare ground under unidirectional wind or 
water stream decays into a chain of barchans after some 
distance of migration (cf. Fig. [5£i,b). Evidence for this 
instability of transverse dunes has been also reported from 
field observations [ill [43[: transverse dunes emerging in 
an area of dense sand cover (e.g. a flat sand beach in a 
coastal area) decay into barchans after migrating a certain 
distance on the bedrock downwind. Therefore, this type of 
dune instability clearly plays a crucial role for the genesis 
and dynamics of barchan dune fields. 

Previous studies of dune genesis have focused on the 
growth of transverse dunes from a sand hill or a flat sand 
bed subjected to a wind of constant direction 1461 1471. l48l 
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I53L |6j. However, the stability of the 
classical transverse dune shape has remained an open issue 
for several decades [l5| . A quantitative study of the long- 
term evolution of a transverse dune was performed recently 
by means of numerical simulations using a model for sand 
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Figure 1: a. Barchan dunes in Westsahara, near 26.51°N, 13.21°W. b. Transverse dunes in Baja California, near 
28.06°N, 114.05°W. In both images, the orientation of the dunes show that the wind is roughly undirectional and blows 
from the top. Images courtesy of Labomar. 
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Figure 2: Evidence for the transverse instability of dunes 
from: (a) numerical simulations 44| — spatio-temporal 
sketch of the profile of a transverse dune at different times, 
which is evolving under constant wind direction. Wind 
blows from the top. The dune is unstable and decays after 
some distance of migration into a chain of barchans; (b) 
water-tank experiments Iff : snapshots of one experimen- 



tal realization of a transverse dune evolving under con- 
stant water stream (flow direction is from the top). An 
initial sand barrier evolves into a transverse dune of width 
~ 3 cm, which then decays into a chain of barchans. Im- 
ages credit: Sylvain Courrech du Pont - Laboratoire MSC 
(Matiere et Systemes Complexes) - UMR 7057 - Universite 
Paris Diderot, Paris, France. 



transport in three dimensions 44[. More precisely, it was 
shown that transverse dunes are unstable with respect to 
any along-axis perturbation in their profile, regardless of 
wind strength, dune size or whether the perturbation is 
random or periodic [441 ]. No matter the type of initial 
perturbation, the transverse dune breaks into a chain of 
barchans, each displaying an average cross-wind width of 
the order of 10 times the original dune height. Indeed, 
the total time of the decay process depends only on the 
dune height and the magnitude of the initial perturbations, 
being independent of the wavelength of the perturbation 

3- 

The question which we want to address in the present 
paper is the following: why does an along-axis perturba- 
tion in the dune profile grow? Here we want to address this 
issue in a quantitative sense performing a linear stability 
analysis of the transverse dune. While numerical simula- 
tions can be used to make quantitative predictions of the 
long-term behaviour of a transverse dune 44j . the linear 
stability analysis refers to the early stage of the perturba- 
tion growth, when the magnitude of the perturbation is 
small [3, EH- We consider an infinitely long transverse 
dune to which a harmonic modulation of small magnitude 
is added. The time evolution of this perturbation is calcu- 
lated by considering a strictly unidirectional wind regime. 

This paper is organized as follows. In the Section 2 we 
give a brief review on the physics of sand transport and 
the formation of transverse dunes. The linear stability 
analysis of the transverse dune is the subject of Section 3. 
Finally, conclusions are presented in Section 4. 

2. The physics of dune formation 

2.1. Formation of transverse dunes from a longitudinal 
sand-wave instability 

A sand plane subjected to a wind of sufficient strength 
is unstable and gives rise to a chain of transverse dunes 

H, S3 El Si, Hi 



56l . l52t l53j . This instability is of hydro- 



dynamic origin. A bump or small hill poses an obstacle to 



2 



the wind, thus producing an upward force that deflects the 
air flow approaching the bump from the upwind. At the 
upper portion of the bump, a downward force keeps the 
flow attached to the surface, thus pushing the flow stream- 
lines closer to each other and enhancing the air shear stress 
on the top of the hill. However, the maximum of the shear 
stress is not exactly at the crest, rather it is shifted upwind 
with respect to the bump's profile 5(1 lj| |[. So, if the flux 
reacted without delay to a variation in wind velocity, then 
maximal erosion would always take place upwind of the 
crest; sand would be then always deposited on the crest, 
and the bump would grow. 

Indeed, there is a transient length needed for the salta- 
tion flux to adapt to a change in wind speed. This length- 
scale is the so-called saturation length [181 Il4|. 
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where Pgrains and Pair are the density of the grains and 
of the air, respectively, and d is the grain diameter. For 
quartz particles (plains = 2650 kg/m 3 ) of average diam- 
eter d = 250 fim [1], L sat for saltation under Earth con- 
ditions (p a i r = 1.225 kg/m 3 ) is approximately 1 m. The 
bump grows only if its width is large enough such that the 
sand flux attains its maximal value upwind of the bump's 
crest. If the hill is smaller than about 20L sa t 0, then it 
is completely eroded and disappears. Conversely, a hill 
that is larger than this minimal size can evolve into a 
transverse dune. The shape of the hill becomes increas- 
ingly asymmetric as the lee slope increases due to depo- 
sition, while the windward side becomes correspondingly 
less steep as a result of erosion. Eventually, the slope at 
the downwind side of the dune becomes so large that the 
flow streamlines there cannot be kept attached to the sur- 
face and flow separation occurs. At the lee, a zone of 
recirculating flow develops which extends downwind up to 
a distance of about 4 — 8 times the dune height [57, 



This so-called "separation bubble" (c.f. Fig. 3) functions 
as a sand trap, since net transport there essentially van- 
ishes [9j. Further, when the lee slope exceeds the angle 
of repose of the sand, 8 C w 34°, the surface downwind 
relaxes through avalanches in the direction of the steep- 
est descent, forming the so-called slip-face, as depicted in 
Fig. [3^,. A sharp brink separates the slip-face from the 
windward side, which has slopes of about 5° — 10° 3^, 31 1. 
Dune migration consists of grains climbing up the wind- 
ward side through saltation, being deposited downwind of 
the crest and thereafter sliding down the slip-face through 
avalanches. 

The transverse dune migrates downwind with a velocity 
v that is roughly inversely proportional to the dune height 
(H). In fact, the larger the size of the dune, the more time 
is needed for a dune to migrate a distance equivalent to its 
width (Lq). The dune turnover time may be understood 
as the time needed for a sand grain, once buried on the 
foot of the slip-face, to reappear on the foot of the dune's 

windward side 58, sum. 



2.2. Barchan dunes 

In the discussion above, we have considered a sand hill 
that is invariant in the direction perpendicular to the wind, 
i.e. an infinitely long transverse sand ridge. If the initial 
surface is not a transverse ridge but rather a small heap 
approximately as long as it is wide, then the longitudi- 
nal instability leads to a barchan dune (Fig. [3j}). Due 
to the scaling v oc 1/H [1], the central slice of the heap, 
which is the one with the largest height, migrates with 
the smallest velocity. In the opposite, the sidemost slices 
advance downwind the fastest, and in this manner a cres- 
cent shape is formed with two limbs pointing in the wind 
direction. Indeed, if the different slices moved uncoupled 
to each other, then they would migrate apart due to their 
different speeds. However, the slices are coupled due to 
lateral sand transport in the direction perpendicular to 
dune migration. This lateral transport occurs on the slip 
face due to gravitational downslope forces arising wherever 
on the slip face the slope exceeds the angle of repose, and 
on the windward side due to the lateral component of the 
wind shear stress [12]. The two factors are negligible near 
the crest where the lateral slope is small and so the central 
slice of a barchan essentially mimicks the cross section of 
a transverse dune of same height. 



3. Linear stability analysis of the transverse dune 

In this Section, we aim to verify the transverse insta- 
bility of dunes through a linear stability analysis. In fact, 
since we do not know the exact analytical form of the 
transverse dune's cross-section in the steady state, i.e. the 
unperturbed solution, it is not possible to perform the full 
stability analysis of the three-dimensional dune. However, 
the problem can be simplified by analyzing the perturba- 
tion at the slip face of the dune, which is where the decisive 
mechanism of lateral transport occurs. 

Following the discussion of the previous section we as- 
sume that the migration velocity of a thin longitudinal slice 
of a transverse dune is essentially inversely proportional to 
its height [l|, or more precisely to the square root of its 
area. So if a perturbation in the area of the slices is ap- 
plied along the transverse direction, then the smaller slices 
move faster, while the larger ones stay behind. In order 
to assure the existence of an instability, it just suffices to 
demonstrate that there is a lateral flux of sand from the 
advanced slices of the dune towards those that are behind. 
If this is the case, then the smaller slices will become even 
smaller, and the larger ones even larger, in such a man- 
ner that their relative speed difference increases, therefore 
enhancing the perturbation. 

We consider an infinite and straight transverse dune, 
which is moving under a constant wind that is strictly 
unidirectional in the x-direction. Let M be the mass of 
this small longitudinal slice of this dune of width dy. So 
M can be approximated as M ss pAdy, where p is the bulk 
density of the dune, taken as a constant in our calculations. 
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Figure 3: (a) Schematic diagram illustrating the longitudinal profile of a transverse dune (full line) and the corresponding 
separation streamline of the wind flow at the lee (dashed line). Inside the so-called "separation bubble", a zone of 
recirculating flow develops, within which net transport vanishes. Both the downwind position x and the dune profile 
h(x) are rescaled by the height of the dune (H). (b) Sketch of a barchan dune of length L and cross-wind width W. 
The central slice of the barchan mimicks the longitudinal profile of a transverse dune of width Lq. 



Since the slip-face can be regarded as a sand trap (c.f. 
Section [5]), the flux leaving a given slice in wind direction 
at the lee side should be negligible. Thus, any change in 
the mass profile occurs due to lateral transport along the y 
axis, i.e. due to the mass exchange between neighbouring 
longitudinal slices. At time t = 0, we add to the area 
of the dune an infinitesimal harmonic perturbation in the 
direction orthogonal to the wind (the y-direction) , such 
that the perturbed area profile can be written as, 



A(y,t) = A [l + eexp(A*)cos(w2/)], 



(2) 



where Aq is the area of the unperturbed transverse dune, 
to and e are the frequency and the initial amplitude of 
the perturbation, respectively, and A is the perturbation's 
growth rate. From Eq. ([2]), it can be seen that the sign of 
A dictates the long-time behaviour of the perturbation. If 
A < 0, then the perturbation decreases and the transverse 
dune is stable, while a positive value of A makes the per- 
turbation increase, which means that the transverse dune 
is unstable. 

The migration velocity of a thin longitudinal slice follows 
the relation v(y,t) = 2c/ y/ A(y, t), where c is a constant 
which depends on the shear stress. Thus, the periodic 
perturbation in the area profile, Eq. (|2"]). induces a pertur- 
bation in the velocity of the form, 



v(y,t) = 



2c 



y/A [l + e exp(At) cos (ujy)] 



Since e is small, we can expand the expression above in 
a Taylor series up to the first order in e, such that v(y,t) 
can be approximated to, 



v(y,t) 



2c 

7Z 



1 — — e exp(At) cos (wy) 



(4) 



Integration over time of Eq. (|4j) gives the position x(y, t) of 
each longitudinal slice, 



x{y,t) = x a 



2c 



-J 



CAq Ay/Ao 



i exp (At) cos (coy) , (5) 



where we assumed, for convenience, that the initial posi- 
tion of the slices is 



x(y,0) = x 



AJA { 



:ecos (toy). 



(6) 



In fact, choosing such an initial profile for x(y,0) means 
adding to the initial transverse dune also a perturbation 
on the xy— plane (besides the perturbation in the height 
given by Eq. @). However, Eqs. (0) and © mean that 
the smaller slices of the initial transverse dune are slightly 
advanced downwind with respect to the larger ones — this 
must occur also if an invariant profile (x(y, 0) = xq) is 
taken instead: due to the scaling v ~ l/H, the smallest 
slices of a transverse dune that is perturbed according to 
Eq. ([5J will move in the front after a certain (infinitesimally 
small) amount of time, such that a perturbation in the 
height also gives rise to a perturbation on the xy-plane 
(Eq. ©). In other words, rather than the specific initial 
profile of x(y, 0), it is the dynamics of lateral transport in 
the course of dune motion which is relevant for the stability 
analysis of the transverse dune. From the profile x(y,t), 
as shown in Fig. 0^, we can obtain the surface equation of 
the slip-face, 



(3) a(x, y) = (x, y, x tan(0 c ) + v t tan (9 C ) - x(y, t) tan (0 C ) 



A^ 



eexp(At) cos (coy) tan (0 C ))> 



(7) 



where vq — 2c/t/A~q. Equation ([7]), which gives the three- 
dimensional profile of the slip-face of the migrating dune 
as a function of time, can be now used in order to compute 
the sand flux along the dune axis. Avalanches occur along 
the slip-face in the direction that makes the smallest angle 
with the vertical direction, i.e. with the downward vector 
—z. The plane of the slip-face is defined by this vector and 
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the normal vector, n = N /\\N\\, where, 
N 



da da 
dx dy 



(1,0,- tan(6» c )) x (0,1, Duj sin {toy) tan (6 C )) 
(tan(6* c ), — Duj sin (toy) tan (6> c ), 1), 



and 



D 



cc 



: exp(A£). 



(8) 



(9) 



In this manner, the vector normal to the slip-face's surface 
is given by the equation, 



W\ 



C 



(sin(0 c ), -Doj sin (ojy) sin {6 C ), cos(6» c )), 

(10) 



where 



C= Jl + D 2 u 2 sin 2 {uy)sin 2 {9 C ). (11) 

So the direction of the flux is given by the normalized 
vector / = where F is determined by a linear 

combination of h and — z, 

~P = -fh~/3z. (12) 

Since F is parallel to the surface of the slip-face, we con- 
clude that F ■ h — 0. Combining this equation with 
Eq. (|T2"j) . we obtain a relation between 7 and j3, 



(3cos(0 c ) 
C ' 



7 — j3z ■ h = ==>■ 7 = f3z ■ h = 
In this manner, Eq. (|12l) becomes, 

^ = ^(sin(0 c )cos(0 c ), 

— Duj sin (ojy) sin (0 C ) cos((9 c ), cos 2 (£? c ) — C 2 



(13) 



(14) 



from which we obtain the normalized vector field of the 
flux of avalanches (see Fig. |4]d). 



/ 



1 



f = 

Duj sin {uiy) cos (9 C ) sin(6* c ), cos 2 (0 c ) — C 2 ). 



:(cos(# c ) sin(6> c 



(15) 



Condition for avalanches — Avalanches will occur in the 
direction of / (Eq. (fT5|) ) wherever the angle 9 defined by 
/ and the xj/-plane exceeds the angle of repose, 9 C = 34°. 
The angle a between / and the vertical is given by the 
equation, 



cos(a) = / • (- 



cos 



-C 2 



^/C 2 ~cos 2 (9 c ) (16) 



cos z 



C 



In fact, the critical angle is measured between the plane 
of the slip-face and the horizontal plane, i.e. 9 — \ — a, 
such that, 



sin(0) = 



VC 2 -cos 2 (0 c ) 



C 



(17) 




Figure 4: (a) Schematic diagram showing the slip-face's 
surface, a(x, y) (defined as in Eq. ([7])) of the transverse 
dune perturbed according to Eq. @. (b) Vector field / 
(Eq. <jT5j) ) of the sand flux along the perturbed transverse 
dune's slip- face (viewed from the top). 



Considering 9 = 9 C + 77, where 77 is of second order in e, we 
can then use the following approximation: 

sin(0) = sin(0 c ) cos(f7) + sin(?y) cos(6> c ) w rjcos(9 c ). (18) 

Calculation of the sand flux — Next we apply the conti- 
nuity equation to calculate the transport of sand through 
avalanches. The mass M(y,t) of an infinitesimally small 
slice of the dune width dy and area A(y, t) evolves in time 
as sand is transported along the axis of the transverse 
dune, i.e. in the direction y, which is perpendicular to 
the surface area of any longitudinal slice. Let us call it 
the velocity vector of the sand flow along the slip-face. 
From mass conservation, we obtain, 

M (y, t + dt) = M (y, t) - {plt(y + ^, t) ■ yA(y, t)dt 



+ plt(y-^-,t)-(-y)A(y,t)dt}. 



(19) 



Since M(y,t) ~ pA(y,t)dy, we can write, 

- P tf(y,t)- d ^ d i)-yA(y,t)}. 



Thus, 



9A(y,t) dlt(y,t) 
P — 57 — = ~P k yA(y,t), 



dt dy 
which leads to the following equation, 

dA(y,t) _ _ du y (y,t) 



c)t 



dy 



(20) 



(21) 



(22) 



where u y is the y-component of the velocity it defining the 
sand flux. The sand flux is proportional to the difference 
between the slope in the y-direction with the critical angle, 
whereas the direction is given by the vector /. By defining 
the velocity vector of the sand flux as it — Br/f, where 
B is the characteristic velocity of the avalanches along the 
slip-face, and using Eq.(15), we obtain, 




-Du sin (uy) cos (9 C ) sin(9 c ) 



Cy/C 2 



cos z 



(23) 



B v/C 2 - cos 2 



cos(6> c ) C 
-Dojs'm(ojy)cos(d c )sm(9 c ) \ 



-B 



C^JC 2 -cos 2 (6» c ) 
Dlu sin {toy) sin(# c ) 
CP ' 



using the values of C and D given by Eqs. (fTTj) and ©, 
respectively, and noting that C ~ + 0(e 2 )), results 



u y(y,t) = -b 

In this manner, 

du y {y,t) = _ B 



ce exp(At)cj sin (coy) sin(# c ) 
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ce exp(At)ui 2 cos (uy) sin(6> c 



(25) 



(26) 



dy ~ A^/Aq 

Combining Eqs. ©, (J22j) and $26$), leads to, 

A Aeexp(At) cos (uy) — A (l + eexp(At) cos (uy))B 
ce exp(At)w 2 cos (uy) sin(9 c ) 



AJA, 



n 



(27) 



Neglecting the terms of second order in e, we obtain, 

. ce exp(At)u 2 cos (uy) sin(# c ) 



j4oAeexp(Ai) cos (uy) = A B- 



AVA^ 



(28) 

and the growth rate of the perturbation can be finally 
written as, 



A = ±0 



S cB sm(9 c ) 
VA^ 



(29) 



Since the positive root in Eq. (|29[) dominates, the pertur- 
bation should always grow in time. Therefore, the linear 
stability analysis corroborates the results from previous 
numerical simulations [44| indicating that transverse 
dunes are unstable with respect to any perturbations 
along their axis. 



4. Concluding remarks 

In conclusion, we have shown, by means of a linear sta- 
bility analysis of the transverse dune, that any perturba- 
tion along the axis of a transverse dune amplifies due to a 
combined effect of the following two main factors: the lat- 
eral transport of sand through avalanches along the dune's 
slip-face, and the scaling of dune migration velocity with 
the inverse of the dune height. It is important to empha- 
size that the growth rate estimated in Eq. (j29|) refers only 
to the initial stage of the dune evolution, where the per- 
turbation can be considered small. As the perturbation 
increases, nonlinear effects become important. Numerical 
simulations showed that, in this regime, the growth rate of 
the perturbations scales with the inverse of dune's turnover 
time, which is the time the dune needs to migrate a dis- 
tance equivalent to its width ji^. It would be interesting 
to conduct a theoretical study of stability where nonlinear 
effects are considered in order to compare the results with 
previously reported numerical simulations (43 |. 

The results of the present paper apply to a transverse 
dune which is moving isolated on the bedrock in an area 
where the wind is strictly unidirectional. Variations in 



G 



wind direction may enhance the transverse instability and 
accelerate the decaying process of a transverse dune into 
barchans. As shown previously, the surface of a large 
barchan may develop small superimposed structures — the 
so-called sand-wave-instabilities — when subjected to 
a secondary wind direction that makes a small angle with 
the primary wind (the one which forms the dune). The sec- 
ondary bedforms developing on the dune's surface, which 
in fact are small dunes with wavelength of the order of the 
minimal dune size, migrate along the dune's crest there- 
after leaving through the limbs. In this manner, large 
amounts of sand can be carried away from large desert 
dunes due to the occurence of secondary winds obliquely 
to the direction of motion. The role of wind trend vari- 
ations for growth rate of the transverse instabilities is an 
open question that remains to be investigated in the fu- 
ture. Furthermore, spatial variations in the influx along 
the direction perpendicular to the wind, e.g. due to the 
presence of dunes upwind, should also contribute to break 
the dune, and must be thus considered in a more detailed 
study of the transverse instability. 

Our calculations elucidate the observation that, when 
the wind blows nearly from the same direction, the dom- 
inant dune shape moving on the bedrock is the barchan, 
and not the transverse dune IB], 44[. In fact, the present 



linear stability analysis shows that real transverse dunes 
moving under unidirectional winds are always unstable, 
independently of the amount of sand on the ground. Our 
calculations do not distinguish between high and low 
sand availability, i.e., the physical mechanisms driving 
the transverse instability exist also when transverse dunes 
evolve in areas of high sand availability, as for example 
on a deep sand bed. Indeed, real transverse dunes are 
never straight, but always display a characteristic sinu- 
soidal shape, which is a fingerprint of the intrinsic insta- 
bility reported here. However, a transverse dune can only 
fully break apart into a chain of barchans if it is moving 
on bedrock in an area with low sand supply }44| . Further- 
more, since the linear stability analysis cannot go beyond 
small deviations from the perfect translational invariance, 
it can not discriminate between a barchanoidal chain, a 
decomposition into individual barchans, or a wavy shaped 
transverse dune with eventually propagating waves. The 
study of the long-term evolution of the transverse instabil- 
ity when the dune is on a sand bed, or when the ground 
is wet or vegetated, might be tackled in the future with 
large scale numerical simulations 44 1. 
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